OFFICE  OF  NAVAL  RESEARCH 
'’.rant  No.  N00014-90-J-1263 

RCT  Project  4133002 - 05 

Technical  Report  #17 


A  SIMPLE  MODEL  OF  ION  SOLVATION  WITH  NON  ADDITIVE  CORES 


by 


W.  R.  Fawcett*  and  L.  Blum** 


Prepared  for  Publication  in  the 
Ukr.  Acad.  (M.  Holovko  festschrif ft) 


AD-A261  080 


mill 


*Department  of  Chemistry,  University  of  California 

Davis,  CA  95616 


and 


mw 
1 


**Department  of  Physics,  University  of  Puerto  Rico 

P.O.  Box  23343 

Rio  Piedras,  P.R.  00931-3343,  USA 


«00 


Reproduction  in  whole  or  in  part  is  permitted  for 
any  purpose  of  the  United  States  Government 


*This  document  has  been  approved  for  public  release 
and  sale;  its  distribution  is  unlimited 


*This  statement  should  also  appear  in  Item  10  of  Document  Control 
Data  -  DD  Form  1473.  Copies  of  form  available  from  congnizant 
contract  administrator. 


U.  REPORT  SECURlT>  CLASSIFICAriON 


2a.  SECUMITY  CLASSIFICATION  AUTHORITY 


2b.  OECLASSIFICATION/OOWNCIUOINO  SCHEDULE 


a.  PERFORMING  ORGANIZATION  REPORT  NUM8ER(S) 

Technical  Report  #17 


REPORT  DOCUMENTATION  PAGE 


10  R£STRlCT;VE  MARKINGS 


60  OFFICE  SYMBOL 
(tf  ip0li<tbl«) 


6a.  NAME  OF  PERFORMING  ORGANIZATION 

Physics  Department 
University  of  Puerto  Rico 


6c  ADDRESS  (Cty.  Star#,  and  ZiPCod*) 

Rio  Piedras 

Puerto  Rico  00931-3343 


8a.  NAME  OF  FUNDING /SPONSORING 

ORGANIZATION  Chemistry 
Office  of  Naval  Research 


8<.  ADDRESS  rcity.  Stitt,  ind  ZIP  Cod*) 

Arlington 

Virginia  22217-5000 


y  y  Title  (/nc/u^  security  CUinh<^tiOfs\ 

A  Simple  Model  for  Ion  Solvation  with 


12  personal  AUTHOR(S) 

W.R.  Fawcett  and  L.  Blum 


8b  OFFICE  SYMBOL 
(tf  ipplicibi*) 

Code  472 


J  O'STRIBUT'ON/ availability  OF  REPORT 


S  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 


7i  NAME  OF  MONITORING  ORGANIZATION 


Tb  ADDRESS  (C/fy,  Srara.  and  ZIP  Cod*) 


9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NlMBER 

RCT  Project  4133002—05 


10  SOURCE  OF  FUNDING  NUMBERS 


PROGRAM  I  PROJECT 

ELEMENT  NO  I  NO 


Non  Additive  Cores 


WORK  UNIT 
ACCESSION  NO 


13a.  type  OF  REPORT 

Paper 


16.  SUPPLEMENTARY  NOTATION 


13b.  TIME  COVERED 
FROM  TO 


14  DATE  OF^REP^RT  {Y**f.  Month.  Diy)  jlS  PAGI^COUNT 


COSATI  CODES 


GROUP  1  SUB-GROUP 


18.  Subject  terms  (Continue  on  reverse  if  necessary  and  idtntify  by  blo<k  numbtr) 


19  ABSTRACT  (Continue  on  r*v*n*  if  necessary  and  idtntify  by  blotk  numb*r) 

The  solvation  i)f  simple  ions  is  studied  with  a  model  with  non  additive  hard 
cores  using  the  ion-dipole  mean  spherical  approximation. 


20.  DISTRIBUTION/ availability  OF  ABSTRACT 
□  UNCLASSIFIEO/UNLIMITEO  □  SAME  AS  RPT  □  O^C  USERS 

Zl  ABSTRACT  security  CLASSIFICATION 

22a  NAME  OF  RESPONSIBLE  NOlViOUAL 

Dr.  Robert  Nowak 

00  FORM  1473, 84  MAR 


SECURITY  classification  OF  TmiS  jaGE 


33  APR  eO'tion  /nay  be  vseo  emausred 
All  other  editions  are  oosoiete 


A  Simple  Model  for  Ion  Solvation  with  Non  Additive  Cores* 


L6 


1 


W.  Ronald  Fawcett 

Department  of  Chemistry 
University  of  California 
Davis,  CA  95616 
U.S.A. 


and 


Lesser  Blum 
Department  of  Physics 
P.O.  Box  23343 
University  of  Puerto  Rico 
RioPiedras,  PR  00931 
U.S.A. 


*In  honor  of  Dr,  Miroslav  Holovko 


DTIC  QUAUTY  IK3PEUTED  3 


AooessioD  For  ..  t 

1  HTIS  GRAJtl 

or 

DTIC  TAB 

□ 

Unannounced 

□ 

Justification — 

By_-_ - 

_ _ ^ 

Distribution/ 


Availability  Ccda^ 
Avail  andyor 
Dlst  Special 


4 


r 


2 


INTRODUCTION 

Ionic  solvation  has  been  the  subject  of  cunent  interest  for  at  least  a  century.^  Modem 
experimental  techniques^  and  the  availability  of  powerful  computers^  have  produced  real  progress 
in  the  understanding  of  this  complex  problem. 

It  is  clear  that  in  many  problems  the  knowledge  of  the  very  complex  details  of  the  close  range 
interactions  are  more  than  what  one  needs.  In  the  interpretation  of  thermodynamic  and  kinetic 
parameters  averages  over  large  domains  are  predominant,  and  therefore,  simple,  yet  controlled 
theories  not  only  provide  a  way  of  fitting  data,  but  also  useful  physical  insight. 

The  success  of  the  Marcus  theory  of  electron  transfer  kinetics^  in  solution  is  a  good  example. 
In  this  theory  the  polarization  of  the  solvent  is  evaluated  from  a  continuum  model  for  the  solvent, 
that  is,  the  structure  of  the  solvent  is  completely  ignored.  The  recent  rather  spectacular  calculations 
of  Chandler  and  his  group^  have  conclusively  shown  that  within  certain  limitations,  this  is  not 
unreasonable  provided  that  the  radius  of  the  activated  complex  is  chosen  in  a  suitable  manner. 

There  are  also  other  examples  where  this  is  true.  It  was  used  by  many  investigators  in  the  past,  but 
in  an  entirely  empirical  way.^ 

In  a  series  of  recent  papers^**®  we  have  discussed  the  use  of  the  Mean  Spherical 
Approximation  (MSA)  in  combination  with  a  very  simple  asymmetric  force  model,  which  has  a 
spherical  repulsive  core.  The  excellent  correlations  for  the  thermodynamic  and  kinetic  parameters 
for  a  large  number  of  solvents  are  more  than  mere  coincidence.  We  believe  that  they  are  due  to  the 
fact  that  the  MSA  is  not  only  a  simple  theory,  but  it  also  has  some  remarkable  variational 
properties,  that  were  discovered  quite  recently.  Liquid  state  theories  like  the  MSA  and  the  HNC 
can  be  derived  as  variational  problems  of  the  free  energy  functional,  which  is  written  in  terms  of 
the  Omstein-Zemike  direct  correlation  function.  This  view  was  introduced  by  Rosenfeld^^‘^3  fw 
systems  of  hard  objects  in  general,  and  hard  spheres  in  particular.  In  that  work  the  relation 
between  the  grand  potential,  the  direct  correlation  function  and  the  scaled  particle  theory  was 
clearly  shown.  This  analysis  was  extended  to  the  case  of  charged  hard  spheres  in  general,  and  to 
the  primitive  model  of  electrolytes  in  particular.*^’^^  it  can  be  shown  that  the  thermodynamics  and 
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the  pair  correlation  functions  can  be  derived  from  a  simple  model  in  which  the  ions  are  replaced  by 
charged  shells  and  the  dipoles  by  charged  dipolar  shells.  We  have  shown  that  the  direct  correlation 
function  in  the  MSA  is  just  the  electrostatic  energy  of  the  shells.  Solving  the  MSA  is  thus 
equivalent  to  an  electrostatics  problem. 

In  the  limit  of  high  density  and  charge  the  Gibbs  energy  and  the  internal  energy  should 
diverge  to  the  same  order  in  the  coupling  parameter  (charge/temperature  ratio)  while  the  entropy 
should  diverge  at  a  slower  rate  and  therefore  the  Gibbs  energy  and  the  internal  energy  coincide  in 
this  limit.  The  mean  spherical  approximation  (MSA)  and  the  hypemetted  chain  approximation 
(NHC)  satisfy  these  bounds,  discovered  by  Onsager  some  fifty  years  ago.  (This  is  not  true  for 
the  Debye  Hiickel  theory,  where  the  entropy  also  diverges).  Therefore,  in  this  limit,  the  HNC  and 
the  MSA  Gibbs  energy,  and  its  derivatives,  coincide.  This  is  a  very  gratifying  feature,  because  the 
HNC,  which  is  the  more  accurate  theory,  is  in  general  difficult  to  solve  numerically,  while  the 
MSA  is  analytical  in  most  cases,  and  of  a  rather  simple  form. 

In  the  asymptotic  limit  the  excess  electrostatic  energy  is  identical  to  the  exact  Onsager  lower 
bound,  which  is  achieved  by  immersing  the  entire  hard  core  system  in  an  infinite  neutral  and 
perfectly  conducting  (liquid  metal)  fluid.  In  the  case  of  the  primitive  model  of  electrolytes,  this  is 
equivalent  to  smearing  the  charge  on  the  surface  of  the  ion.  It  was  also  shown^'***®  that  the  direct 
correlation  function  is  equal  lo  the  interaction  potential  of  two  "smeared"  ions  of  radius  Ri,  Rj 
where  Rj  is  the  radius  of  ion  i  plus  some  screening  length. 

The  Onsager  process  of  introducing  the  infinite  conductor,  naturally  decouples  all  the 
different  components  in  the  system  which  may  differ  in  size,  shape,  charge  distribution;  and  relative 
orientation  in  space.  As  a  result,  the  variational  free  energy  functional  in  the  high  coupling  limit 
diagonalizes,  and  the  mathematical  solution  of  the  asymptotic  problem  is  given  in  term  of  the 
geometrical  properties  of  the  individual  particles  in  the  system. 

The  MSA  Gibbs  energy  provides  an  appropriate  functional  from  which  the  charge  part  of  the 
direct  correlation  function  is  obtained  by  functional  differentiation. 


There  is  an  actual  numerical  test  case  that  proves  the  electrostatic  analogy  of  the  MSA.  The 
very  high  charge  and  density  limit  for  the  one  component  plasma  (point  ions  in  a  neutralizing 
background)  was  investigated  using  the  HNC  equation  by  Ngl^,  in  a  tour  de  fourcer  numerical 
calculation.  The  same  results  were  obtained  by  Rosenfeld  et  al.  using  the  analytical  solution  of  the 
MSA.  18 

A  second  point  is  the  use  of  spherical  cores  to  represent  non  spherical  objects.  This  also 
depends  on  the  properties  of  interest.  For  the  equation  of  state,  it  works  even  for  very  non 
spherical  objects.  There  are  numerous  examples  of  sphericalization  in  the  literature  of  molecular 
fluids.  A  notable  example  is  given  in  the  recent  work  of  Williams  et  al.l^  The  median  diameter  of 
a  non  spherical  object  yields  a  surprisingly  good  equation  of  state,  not  only  for  hard  dumbbells  but 
also  for  real  gases  at  extremely  high  pressures.^O*^!  Other  examples  are  found  in  the  extensive 
studies  of  Nezbeda  and  his  group22.23  on  the  structure  of  molecular  liquids. 

The  analytical  solution  of  the  MSA  for  mixtures  of  ions  and  dipoles  with  a  non  spherical 
sticky  potential  was  obtained  some  time  ago2^*26  for  the  case  of  a  sphere  with  a  point  dipole  and  a 
sticky  potential  of  dipolar  symmetry,  which  represents  the  specific  non  electrostatic  interactions. 

In  this  woric  it  was  shown  that  the  contribution  of  the  sticky  potential  was  merely  a  shift  in  the 
polarization  parameter,  X.  This  work  has  been  continued  to  include  short  ranged  potentials  of 
higher  (tetrahedral)  symmetry.22 

In  this  spirit  we  have  recently  analyzed  the  solvation  thermodynamics  of  monoatomic 
monovalent  ions  in  water  and  in  a  number  of  polar  solvents.2.l0>28  The  basic  idea  comes  from  the 
observation  that  to  a  first  approximation  the  corrections  due  to  the  molecular  nature  of  the  solvent 
result  in  a  change  in  the  screening  length  from  lA*  to  l/r  +  Xoi/r  0$,  the  shift  being  different  for 
each  ion,  according  to  its  size;  because  of  geometrical  packing  considerations  the  number  of 
solvent  molecules  that  can  be  the  first  neighbors  of  a  given  ion,  i,  vary  according  to  the  ratio  Oi/Os. 

Consider  now  a  realistic  model  of  the  solvent,  with  a  soft  repulsive  core,  and  of  non-spherical 
shape.  Then,  the  distance  of  closest  approach  of  a  certain  ion  should  be  very  different  for  the 
anion-solvent  pair  and  the  cation-solvent  pair,  even  if  the  unsolvated  anion  and  cation  are  of  the 


same  size.  The  correct  statical  mechanical  model  for  such  a  situation  is  the  mixture  with  non¬ 
additive  cores  discussed  in  the  literature.^^  A  simple  way  of  accounting  for  this  is  to  assume  that  at 
low  concentration  the  effective  radius  of  the  solvent  is 


^si  -  fs  ^si 


(1) 


The  Model 

In  the  MSA  for  a  system  of  spherical  ions  and  dipoles  of  arbitrary  size^®,  the  Gibbs  energy 
of  solvation  in  the  limit  of  infinite  dilution  is 


,  .  NaUiSflll  A  .  i'l  _i 

Sic  Eo  V  es 


)  ri+50 


(2) 


Here,  z,  is  the  valence  of  the  ion,  e© ,  the  fundamental  charge,  Es,  the  static  dielectric  constant  of 
the  solvent,  Eq.  the  permittivity  of  free  space,  rj,  the  radius  of  the  ion,  and  Nq,  the  Avogadro 
constant.  5°  is  equal  to  rj/X,  where  rs  is  the  radius  of  the  solvent  (the  effective  radius)  and  the 
MSA  polarization  parameter,  X  is  calculated  from  the  dielectric  constant  of  the  pure  solvent  £5 
using  the  Wertheim  relationship*** 

XMX+1)4  =  16  Es  (3) 


Although  this  value  of  the  polarization  parameter  may  be  appropriate  for  estimating  dipole-dipole 
interaction  energy,  it  gives  poor  estimates  of  the  Gibbs  solvation  energy  when  this  is  assumed  to 
equal  Gid.^’^®’^^  Therefore,  and  in  view  of  the  above  discussion  we  change  eq.  (2)  to 


where 


Gid(C+) 


NaijiCcH  .  U  _±_ 

SjtEo  V  csj  (n+Sc) 


5c 


rs+Asc 


and  for  cations,  and  similarly,  for  monovalent  monoatomic  anions 
Gi,(A-)  =  -  A  .  L\  _J_ 

SttEo  1,  EsJ  (H+Sa) 


(4) 

(5) 

(6) 
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where 

5a  =  (7) 

This  separation  and  parameterization  produces  a  much  more  sensible  and  physical  model  of  Gjd  for 
ions  in  a  number  of  solvents. 

The  second  contribution  to  the  Gibbs  solution  energy  comes  from  repulsive  dipole-dipole 
interactions^S.*  and  is  given  by  the  equation 


Gdd  = 


(zieo)^(es-l)^ 
647t  eoEs  (ri  +5i)2 


4ri  +  6 

lTxT7 


2(X+1)2  J 


(8) 


This  expression  may  be  simplified  considerably  when  one  considers  the  range  of  values  typical  for 
X.  For  water  whose  dielectric  constant  is  78.3  at  25  ®C,  X  is  equal  to  2.65.  The  corresponding 
value  of  the  ratio  X(X+3)/2(X+l)2  which  appears  in  the  denominator  of  the  term  in  square  brackets 
in  eq.  (8)  is  0.56.  This  is  negligible  with  respect  to  the  value  of  £$  and  may  be  neglected.  The 
ratio  appearing  in  the  numerator  of  this  term,  namely,  (3X,+2)/(X+l)  is  equal  to  2.7.  Thus,  the 
expression  for  Gdd  niay  be  rewritten  as 


647teo  V  Esj  (ri+5i)2 


(9) 


The  simplified  expression  has  the  advantage  that  the  parameter  X  does  not  appear  on  its  own  but 


only  in  5i.  This  fact  is  important  in  applying  the  expression  to  experimental  data.  Finally,  on  the 
basis  of  previous  woric^S,  only  a  fraction  9  of  the  dipole-dipole  term  is  used  in  estimating  the 


Gibbs  energy  of  solvation  so  that  the  final  expression  for  Gs  becomes 

Gs  =  Gid  +  0  Gdd 


(10) 


The  expressions  for  the  corresponding  entropies  are  obtained  from  the  temperature 
derivatives  of  the  above  Gibbs  energies.  For  the  ion-dipole  term,  the  result  for  a  cation  is 
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1  de. 


dSc' 

(  Es) 

dT 

(11) 


StcCo  Les^  (ri+5c)  (ri+5c 
where  d5c/dT  is  the  temperature  derivative  of  8c  which  is  obtained  fix)m  the  temperature  derivative 
of  Xc-  Comparing  eqs.  (4)  and  (1 1),  it  is  easily  shown  that 


Sid(C+)j'l  -  = 


Gid  (C-*-)  ^  8n  EoGid^fC-*-)  dSc 
e,2  dT  -  No(zieo)2  dT 


The  corresponding  expression  for  an  anion  is 


(12) 


Sid(A-)^l 


Gjd  (A~)  dCs  Stc  EoGid^(A‘)  dS/y 
dT  *  No(zieo)2  dT 


(13) 


From  these  equations  it  follows  that  two  additional  parameters  are  required  to  estimate  the  entropy 
of  solvation  of  a  salt,  namely,  the  temperature  derivatives  of  the  polarization  parameters  for  the 
cation  and  anion. 

From  the  expression  for  Gdd.  one  obtains  the  following  equation  for  the  corresponding 
entropy: 


c. 2  dCs//,  l\.2.75idX/^_,  „„^.^  25idX/,  .c>v 

Sdd  —  Gdd  2  dx  '  *  1  \  dT  *  ■s  dT 

The  derivative  dX/dT  may  be  found  from  the  Wertheim  equation  from  which  one  obtains 


dX 

dT 


X(l+X)]  _1_^ 
L  3X+1  J  2es 


(15) 


Assuming  that  the  fraction  0  is  independent  of  temperature,  the  resulting  expression  for  the 
entropy  of  solvation  is 

Ss  =  Sid  +  0  Sdd  (16) 

This  model  is  now  examined  with  respect  to  data  for  monoatomic  nnonovalent  ion  solvation  in 


(14) 


water. 


Results  and  Discussion 


In  order  to  apply  the  above  model  to  experimental  data  one  must  choose  radii  for  the  ions  and 
an  extrathermodynamic  assumption  for  extracting  single  ion  solvation  energies  from  experimental 
values  for  the  salts,  namely,  the  alkali  metal  halides.  Two  sets  of  ionic  radii  give  good  fits  between 

experiment  and  theoryl®,  namely,  the  values  given  by  Pauling*  and  those  extracted  from  neutron 

3S 

and  X-ray  diffraction  experiments.*  The  Pauling  values  are  used  here  simply  because  they  are 
more  familiar.  The  extrathermodynamic  assumptions  used  to  separate  experimental  quantities  into 
those  for  cations  and  anions  were  examined  in  detail  for  aqueous  systems  by  Conway.*  Values 
of  the  Gibbs  energy  and  entropy  of  solvation  for  the  alkali  metal  cations  and  halide  anions 
estimated  by  Conway  are  summarized  in  Table  1  together  with  the  Pauling  radii  for  these  ions. 

In  our  previous  analysis  of  the  Gibbs  solvation  energy^>^0,  the  contribution  of  dipole-dipole 
interactions  was  ignored  ( 0  =  0  in  eq.  (10)).  Analysis  of  the  data  for  the  alkali  metal  halides 
showed  that  the  best  values  of  5c  and  5a  (eqs.  (4)  and  (6))  assuming  Pauling  radii  for  the  ions 
were  82.4  and  18.0  pm,  respectively.^^  The  quality  of  the  ftt  between  theory  and  experiment  is 
really  excellent,  the  standard  deviation  for  the  calculated  values  being  1.5  kJ  mol*^  that  is, 
significantly  less  than  1  percent.  This  result  is  illustrated  in  Fig.  1  in  terms  of  a  plot  of  -Gs  against 
l/(ri+5i).  It  is  apparent  that  the  values  of  Gs  fit  very  well  on  a  straight  line  in  the  region  where  data 
exist  but  that  the  slope  of  the  best  straight  line  through  these  data  is  slightly  higher  (75.2  kJ  nm 
mol'^)  than  that  predicted  by  theory  (69.45  kJ  nm  mol"').  As  a  result,  an  ion  with  infinite  radius  is 
predicted  to  have  a  Gibbs  solvation  energy  of  35.7  kJ  mol*‘  when  this  simple  model  is  applied. 
This  can  be  attributed  to  the  failure  to  account  for  dipole-dipole  interactions  which  constitute  a 
positive  contribution  to  Gs. 

The  same  data  were  re-examined  on  the  basis  of  eq.  (10)  using  the  expression  for  Gdd 
given  by  eq.  (9).  The  fit  was  carried  out  so  that  the  resulting  estimates  of  Gij  and  Gdj  would  go  to 
zero  for  an  ion  of  infinite  radius.  The  resulting  best  values  of  6t  and  8\  were  much  smaller  than 
those  obtained  by  the  simpler  analysis,  and  are  recorded  in  Table  2.  The  corresponding  estimates 
of  Gid  are  plotted  against  l/(ri+6i)  in  Figure  2.  As  expected  these  estimates  are  larger  in  magnitude 
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than  the  experimental  values  of  Gs  also  plotted  on  the  same  graph.  Both  plots  are  linear  in  l/{ri+5i) 
and  have  zero  intercepts  within  the  experimental  standard  deviation.  Estimates  of  0G(jd  for  the 
same  nine  ions  are  also  plotted  in  this  figure.  For  the  cations,  this  quantity  varies  from  a  high  of 
41.2  kJ  mol-*  for  Li+  to  a  low  of  22.7  kJ  mol**  for  Cs+,  and  is  approximately  eight  percent  of  the 
magnitude  of  Gid-  In  the  case  of  the  anions  it  represents  a  slightly  higher  fraction  of  the  total  Gibbs 
energy  and  varies  from  a  high  of  42.7  kJ  mol**  for  F*  to  a  low  of  30.8  kJ  mol**  for  I*  ion.  The 
values  of  0Gdd  are  also  approximately  linear  in  l/(ri+5)  with  an  intercept  corresponding  to  zero  for 
an  infinitely  large  ion.  The  value  of  5  used  in  estimating  Gdd  is  that  estimated  on  the  basis  of  the 
Wertheim  equation,  namely,  51.6  pm.  When  the  contributions,  Gjd  and  0Gdd  are  added  one 
obtains  estimates  of  Gs  which  agree  with  those  obtained  from  experimental  data  to  within  a  few  kJ 
mol-*.  Values  of  the  parameters  used  in  our  calculations  including  the  fraction  0  are  summarized  in 
Table  2. 

A  successful  model  for  ionic  solvation  must  also  be  able  to  estimate  correctly  the  entropy  of 
solvation.  In  order  to  make  this  calculation  one  must  have  estimates  of  the  temperature  derivatives 
of  the  parameters  6c,  5a  and  X,  The  latter  was  calculated  using  eq.(15)  which  gives  a  value  of 
-2.48  X  10-*  K-*  for  dVdT  on  the  basis  of  the  dielectric  properties  of  water.  Assuming  rj,  the 
radius  of  a  water  molecule,  is  equal  to  137  pm  the  corresponding  value  of  d5/dT  is  0.048  pm  K**. 
Values  of  5c  and  5a  were  determined  from  the  single  ion  estimates  of  Ss  given  by  Conway.  *3  In 
the  case  of  cations  on  the  basis  of  eqs.  (12)  and  (16),  one  may  write 


where 


and 


Y  = 


X  = 


Sff  eoGid2(C-»-) 
No(zieo)2 


Y  =  Ss(C^)  + 


Gid(C-^)  ^ 

es2  dT 


ss 


(17) 


(18) 


(19) 


Sc(C+)  being  the  solvation  entropy  of  the  cation.  Similar  equations  can  be  written  for  anions  using 
eqs.  (13)  and  (16).  A  plot  of  Y  against  X  for  both  ions  is  shown  in  Fig.  3.  These  data  were  fit  to 
a  straight  line  passing  through  zero  (one  parameter  fit)  as  shown,  the  resulting  slope  which  equals 
both  d5s/dT  and  dS^dT  being  0.032  pm  K'^.  It  is  important  to  note  that  the  present  analysis 
suggests  that  these  temperature  coefficients  are  equal  in  contrast  to  our  previous  conclusion  based 
on  an  analysis  which  ignored  dipole-dipole  interactions.  The  estimated  error  in  the  slope  is  not 
large  being  about  4%.  It  should  also  be  noted  that  the  magnitude  of  d5c/dT  and  d5A/dT  is  about  65 
percent  of  that  found  on  the  basis  of  the  Wertheim  equation  (eq.  (15)). 

Using  the  above  coefficients  values  of  Sid  were  estimated  for  the  alkali  r.ijtal  and  halide  ions 
and  are  plotted  in  Fig.  4  together  with  estimates  of  Sdd  and  experimental  values  of  Sg  against  the 
experimental  Gibbs  solvation  energy  Gs.  It  is  clear  that  the  entropy  of  solvation  is  approximately  a 
quadratic  function  of  the  Gibbs  energy  as  suggested  by  the  MSA  model. 


a 
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